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Abstract: In order to understand thermodynamical properties of A^ D-branes with 
>- ■ chemical potentials associated with i?-syninietry charges, we study a one dimensional 



large N gauge theory (bosonic BFSS type model) as a first step. This model is 
obtained through a dimensional reduction of a 1 + D dimensional SU{N) Yang- Mills 
theory and we use a 1/D expansion to investigate the phase structure. We find three 
(^ ■ phases in the ft — T plane. We also show that all the adjoint scalars condense at 

Q ■ large D and obtain a mass dynamically. This dynamical mass protects our model 

from the usual perturbative instability of massless scalars in a non-zero chemical 
potential. We find that the system is at least meta-stable for arbitrary large values 
of the chemical potentials in D — )■ oo limit. We also explore the existence of similar 
condensation in higher dimensional gauge theories in a high temperature limit. In 
2 and 3 dimensions, the condensation always happens as in one dimensional case. 
On the other hand, if the dimension is higher than 4, there is a critical chemical 
potential and the condensation happens only if the chemical potentials are below it. 

Keywords: 1/N Expansion, Confinement, Brane Dynamics in Gauge Theories, 
M(atrix) Theories. 
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1. Introduction and summary 

Supersymmetric large N gauge theories with 16 supercharges, which describe A^ D 
brane dynamics, are fundamental theories in the study of the string theory. How- 
ever the thermodynamical properties of these theories with /^-symmetry chemical 
potentials have not been studied sufficiently in terms of the gauge theory. In the 



strong coupling regime, the gravity analysis is valid through the gauge/gravity cor- 
respondence, where the i?-symnietry charges in the gauge theories will correspond 
to electric charges or angular momenta of the dual gravity systems [1, 2, 3]. Impor- 
tantly, it has been shown that gravity systems with these charges have rich phase 
structures. For example, in asymptotically AdS space, RN black hole/hairy black 
hole type transitions were found [4, 5, 6, 7] and their applications to condensed mat- 
ter physics have been investigated [8, 9, 10]. Besides, in higher dimensional black 
holes with angular momenta, Meyers-Perry black hole/black ring type transitions 
happen [11]. Therefore we can expect that similar phase structures also appear in 
the supersymmetric gauge theories with corresponding chemical potentials. 

However the gauge theories have not been studied sufficiently, since the presence 
of the chemical potentials makes the analysis difficult even in a weak coupling. The 
reason is as follows. Fields in the gauge theories are typically massless except in some 
special situations, e.g. A/" = 4 SYM on i? x S*^. (Several groups studied the chemical 
potentials in this case [12, 13, 14, 15, 16, 17, 18, 19, 20, 21].) Then the massless 
scalars coupled to the chemical potentials will be unstable around the trivial vacuum 
and perturbative techniques do not work. If these scalars are free, this instability is 
inevitable and the theory is destabilized by the chemical potentials^. On the other 
hand, interactions involving the scalars may remove this instability. For example, 
in a massless g\(t)\'^ theory, if we turn on a finite chemical potential /i for a f/(l) 
rotation, the potential becomes — /i^|0|^ + (7|0|^. Although it causes the shift of the 
vacuum, the theory is still stable ii g > Q. As in this example, the understanding 
of the interactions in the gauge theories is essential to investigate the finite chemical 
potentials. 

In this article, in order to consider this problem, we study the following one 
dimensional large N gauge theory (a bosonic BFSS [23] type model). 






S= f dtTi\ > ■;fDnrM^-> ' -\Y\Y']\Y\Y'] \. {l.V 



Here Y^ are SU{N) adjoint scalars. The covariant derivative is defined by DqY^ 



^A naive regularization in free super Yang-Mills theory through an analytic continuation of the 
chemical potential to complex was proposed in [2]. In this article, we use a different analysis and 
do not consider this regularization. The imaginary chemical potential is considered in lattice gauge 
theory also to avoid the fermion problem in QCD [22]. 



dtY^ — i[Ao, Y^] and Aq is an SU{N) gauge field. This model can be regarded as a 
dimensional reduction of a 1 + D dimensional pure Yang-Mills theory. The model is 
invariant under the SO{D) rotation of F^, which is related to the /^-symmetry in the 
supersymmetric gauge theory. We will consider chemical potentials associated with 
f/(l)L^/2j (- S0{D).'^ Note that the charges of these f/(l)L^/2j correspond to the 
transverse angular momenta on \_D /2\ planes in the context of the rotating D-branes 
[2, 3, 24, 25, 26, 27, 28, 29, 30]. 

We can obtain the model (1.1) with Z) = 9 from N Dl-branes on a spatial circle 
as follows [31]. We can choose the boundary conditions of the fermions for this circle 
as anti-periodic. Then all the fermions obtain a mass, which is proportional to the 
inverse radius of the circle. Thus, if the radius is sufficiently small, we can ignore 
all the fermions, KK-modes and long string states and the theory is reduced to one 
dimension. Then the SO{D — 1) symmetry is enhanced to SO{D) and we obtain 
(1.1). Note that it is known that the gauge/gravity correspondence is not valid in 
such a small radius, since the effective coupling becomes small [31]. Thus we cannot 
compare our results with the gravity directly. Our results should be regarded as a 
weak coupling continuation of the gravity results. (However, as far as I know, there 
is no known gravity result of the thermodynamics of this system with the chemical 
potentials.) 

The thermodynamics of the model (1.1) without the chemical potential has been 
investigated by using a 1/ D expansion [32]. Especially it has been revealed that a 
non-trivial stable vacuum exists in which the adjoint scalars condense and obtain a 
mass. We generalize this analysis to the finite chemical potentials. We will show 
that, in the D — )> oo limit, the (meta-)stable condensed vacuum always exists for 
arbitrary values of the chemical potentials. This result is similar to the g\(t)\^ theory 
and we can conclude that the commutator-squared interaction stabilizes our model 
in the large D limit. However, in a large but finite D case, the 1/ D expansion does 
not converge for large chemical potentials and it is unclear whether the system is 
stable or not there. 

We investigate the phase structure also and show the existence of three phases 
in the ^ — T plane, which is summarized in Figure 1. 

In a single chemical potential case, we find several saddle points in a high temper- 



^Here [• • • J denotes the floor function which maps a number to the integer part. 



ature regime. In this regime, since only the zero-modes of the Matsubara frequencies 
are dominant, the model (1.1) will reduce to an effective zero-dimensional matrix 
model (bosonic IKKT type model [33]). In this model, the chemical potential in- 
duces a CS like interaction. (However, its coefficient is real and different from the 
ordinary CS term in [34].) Therefore, (imaginary) fuzzy sphere like solutions exist 
as complex saddle points. However, physical interpretation of these solutions is yet 
to be explored. 

We also explore the existence of the condensed vacuum in the d dimensional large 
A^ gauge theory to understand the stability issue. '^ Since the analysis is difficult 
in general, we restrict our study to a high temperature regime and consider the 
D — 7- oo limit. Then we will see that the existence of the condensation depends on 
the dimension d. In rf = 2, 3 cases, the condensation always happens for arbitrary 
values of the chemical potentials as in the d = 1 case. On the other hand, in li > 4 
case, a critical chemical potential /Xc exists and the condensation happens only if all 
the chemical potentials are below ^f.. These results are summarized in Figure 4. 

The same condensations will happen in the supersymmetric gauge theories also, 
since we can ignore the fermions in the high temperature limit. Thus our results may 
be related to the analysis in the strong coupling regime through the gravitational 
study of the rotating D brane geometries [3] . However the application of our large 
D results to the supersymmetric gauge theories are not robust particularly in a high 
chemical potential regime and further analyses are necessary. 

The organization of this article is as follows. In section 2, we introduce the 
chemical potentials to our model (1.1) and, by using the large D expansion, we derive 
an effective action. We will analyze this effective action in two regimes. One is a 
low temperature and low chemical potential regime. Another is a high temperature 
regime and/or a high chemical potential regime. We will study the phase structure 
of the first regime in section 3. In section 4, we investigate the nature of the second 
regime. We will also discuss the fuzzy sphere like saddle points in this section. In 



■'The 1/ D expansion analysis of d dimensional gauge theory without the chemical potentials for 
d = {) (bosonic IKKT model) is considered in [35] and for d = 2 (two-dimensional gauge theory on 
T^) is considered in an ongoing work [36]. 



section 5, we explore the existence of the condensation in higher dimensional gauge 
theories. In section 6, we conclude with a discussion. 

2. One-dimensional gauge theory with chemical potential 
2.1 Chemical potential 

We consider the large A^ gauge theory (1.1) with chemical potentials associated with 
the global f/(l)L^/^-l transformations. It is convenient to rename the D adjoint scalar 
F^as 

X'^Y' {I = l,...,[D/2\), 

W' = Y' {1= [D/2\+l,...,2[D/2\), 

Y = Y (If D is even, this one does not exist.), (2.1) 

and define complex scalars as $^ = {X^ + iW^)/^/2, {I = 1, . . . , [D/2\). Then the 
action (1.1) is invariant under the global f/(l)l-'^/^J transformation: ^^ — )> e^^'^^. In 
the context of the rotating D-brane, this charge corresponds to the angular momen- 
tum on the {X^ , W^) plane. 

Now we introduce the chemical potential fij for these U{1) transformations. 
It has been studied in [13, 37] that, in the Euclidean path integral, the chemical 
potential causes the modification of the action as 

D,^ Do -1^1 {I = l,...,[D/2\), (2.2) 

in the kinetic term of $^. Without loss of generality, we can take fij positive. This 
modification gives rise to negative mass term — /i^|$^p and the system might be 
destabilized. In addition, the usual perturbative analysis does not work. 

In the following sections, we will investigate the thermodynamics of the gauge 
theory by using a 1/D expansion [32, 35], in which the number of the adjoint scalars 
is regarded as large. In order to keep the contribution of the chemical potentials 
finite under the large D limit, first we consider a simple situation. 



jj,i = H (/= I,--- ,D), 

/x. = {t = D + l,--- ,[D/2\), (2.3) 



and D scales linearly with D. If we consider general chemical potential /i/, each con- 
tribution will appear as sub-leading in the 1/ D expansion. Thus in order to evaluate 
them consistently, we have to derive other sub-leading terms also and compare with 
them. We will study it in section 3.3 and section 4.3. 

The validity of the application of the 1/ D expansion to Dl-brane case {D = 9) 
is not a priori obvious, since D is large but finite. However in the case of the 
zero-chemical potential, the 1/D expansion [32] reproduces the numerical results 
in [31, 38, 39, 40, 41] quantitatively. In fact the leading large D results even for 
D = 2,3 reproduce the essential qualitative properties. Thus we expect that the 
1/D expansion is also meaningful in the model involving the chemical potentials. 

2.2 1/D expansion and computation of effective action 

In this subsection, we will integrate out the adjoint scalars in the action (1.1) with the 
chemical potential (2.2) and derive an effective action by using the 1/D expansion. 
We evaluate only the leading order of this expansion in this section. We will discuss 
the contribution from the next order in section 3.2 and 4.2. In section 3.2, it will 
turn out that inclusion of these corrections does not change the nature of the phase 
structure in a low chemical potential regime. On the other hand, in section 4.2, we 
will show that the 1/D expansion does not converge in a very high chemical potential 
regime. 

According to [32], we rewrite the path integral by employing an auxiliary field 
Bab as"^ 



Z =M f VBVAoVY'V<l?^e^^^^'^'"^^'^\ 



S{B,Ao,Y,<l>)= I dt 



:BabMj,Bcd 



4^2 
D D 

2 



+ Y,^\^HDo-A + Bab)H+ J2 lY:{-Dl, + Bab)Y^ 

-f=l i=2D+l 

(2.4) 



Again we have used the notation Y^ for the adjoint scalars, which do not couple to 
the non-zero chemical potentials. Here 1/Af = fVB cxp (^j dtBabM^^^^Bcd/{'lg'^)) is 



*The definition of Bab is different from [32] by an imaginary factor "z" . 



a numerical factor. M^^^^^ is the inverse of Mab,cd, which is defined by 

Mab,cd = --|Tr[A„, AJ[Afe, AJ + (a ^ 6) + (c ^ rf) + (a ^ 6, c ^ rf)}. (2.5) 

Here A^ (a = 1 . . . A^^ — 1) is a generator of SU{N). They satisfy Mab,cdM~J^^f = 
{Sae^bf + SafSbe)/'^. Properties of the matrix Mab,cd is summarized in appendix A 
of [32]. We can reproduce the original action (1.1) from (2.4) by substituting the 
solution of the classical equation of motion for Bab'- 



\ 1=1 i='}n^i I 



(2.6) 



Now we can formally integrate out $^ and F*, since the action (2.4) is quadratic 
in them. Then we obtain an effective action for Bab and Aq, 



Seff{B,AQ 



dt 



Ag 



-BabMj^^^Bcd 



k D 
+ ^ logdet {-Dla, + Bab) + ^ logdet (-(Do - A + Bab) , 



kiD 
T 



(2.7) 



where we have defined normalized ratios: ki = {D — 2D)/D and k2 = 2D/D, which 
satisfy ki + k2 = 1. As in [32], we investigate this model by taking the large D and 
large A^ limit such that D -^ oo, N -^ oo and fif — > with fixed A = g'^ND = XD. 
(In our case, we also take D — )■ oo with fixed k2 = 2D/D). Then the first term and 
the logdet terms in (2.7) will be comparable in this limit. As a result, this model 
will have a non-trivial saddle point Bab = '^o^ab, which we will confirm later^. Here 
Ao is a time independent constant. From equation (2.6), indeed this saddle point 
gives a condensation of the adjoint scalars 



D 



D 



N 



2E(Tr$'^<^')+ E (Trrr) = ^A2, 
where we have used a relation M~^^^5cd = j^^ab [32]. 



(2.^ 



^The large N limit is not necessary to derive this saddle point. Actually, even at finite N, we 
can calculate some physical quantities by taking D —>■ oo and g'^N — >■ with fixed g^ND [32, 35]. 
However our interest is in the large N limit of (1.1) and we do not consider finite N effects in this 
article. 



To proceed, we write the Bab as the sum of a constant trace piece and the rest 
as 

Bab{t) = A^6ab + gbab{t), (2.9) 

where bab{t) satisfies J dt baa{t) = 0. Note that we can ignore the interactions between 
bab and $^ and F* in the leading order of the 1/D expansion [32]. The contributions 
of the interactions will appear in the next order as we will see in section 3.2 and 4.2. 
Here we consider a gauge fixing of Aq. It is convenient to take the constant 
diagonal gauge: Aoij(t) = aidij. As a result the effective action will be described by 
gauge invariant Wilson loop operators, 

1 1 ^ 

Un = — Tre^"^''^"'^* = — ' 
" N N 

1=1 

This gauge fixing gives rise to a Faddeev Popov determinant [42] 



1 1 

— Tre*" -^o" ^i"^* = _ \^ e'"'^"' . ('2.10) 



VA 







J]rfa,e-^^^, Spp = iV2 5^ i|M„|2. (2.11 



Now we integrate out bab and obtain the effective action for the condensation A 
and the gauge field {«„}. 

ge//(A,K}) _ /3A4 1 y.|tX„p 

n=l 

+ ^ [A;iTrlog {-DI + A') + A;2Tr log (- {Do - fif + A')] . 

(2.12) 

Here the first term is a classical term from the first term in (2.7). The second term 
comes from (2.11) and indeed it is 1/D order. We keep it here, since this term will 
be dominant in a low temperature regime and more significant than other 0{1/D) 
terms. The integral of bab gives a numerical factor which cancels out A/" in (2.4) 
except the contribution of the integral of the constant trace peace A. 

Let us evaluate the terms in the second line of (2.12). In the momentum space, 
the quadratic term of <l>^ in (2.4) can be written as 



f 2nn _ 



9 V mj " mjj 




(2.13) 



Then the eigenvalues of this matrix are calculated as 

l-nn 



/3 



[aj — q;j) ± iyU J + A . 



(2.14) 



By using this result, we calculate [13, 42] 
Trlog(-(Do-^)' + A2) 

4e 



n,t,j 



log 



2nn 



- {aj -ai)+i^] + AM + lo, 



27rn 



- (ttj - ai) -ill] + A 



:5^1og[Are^^(l 



-/3A+J/3(ai-a^-)-/3M 



)(i 



~l3A-ip{ai-aj)+PiJ.'\ 



n-,i,J 



=Ar2 log;^ + N^/3A -N^y^ -e~"^^ (e"^^ + e""'^'^) 



I ""«! ) 



(2.15) 



n=l 



where we have ignored O (l/N^) corrections. A/" is a irrelevant constant factor and 
we will ignore it from now. Similarly we evaluate the term from Y^ integral also. 
Then the effective action (2.12) becomes 



oo I 19 

7/ 



%/(A,K})_ /3A^ I 1 y- 

8A D^ n 



DN^ 



n=l 



^^-t 



n=l 



,-n/3A 



h + k2 



^-nl3fi _|_ ^nl3fi 



Ur. 



n 



(2.16) 



Note that, during this derivation, we have assumed a relation A > /i. If this in- 
equality is not satisfied, tachyonic modes will appear and the path integral is not 
well defined^. Later we will show that this assumption is ensured. 

From now we evaluate the effective action (2.16) and investigate the phase struc- 
ture and the condensation of the scalars. It is convenient to integrate out A first 
and derive the effective action for the Wilson loops. In order to do it, we analyze 
the saddle point equation for A^, ^ 



A3 1 °° 

2A 2 ^ 



h + k2 



-nPn _j_ gn/3^ 



\Ur. 



0. 



(2.17) 



^ Is — fi case is also subtle, since massless modes will appear and the effective action will be 
non-local. Indeed it will happen at T = 0. We will come back to this problem later. 

^When we derive the saddle point equation (2.17), we deviate the effective action with respect 
to A^ (not A), since A^ is the correct variable as in (2.9). 



We will solve this equation in two different regimes characterized by the values of 
the Wilson loop operators. One is the low temperature and low chemical potential 
regime (l-u^l ~ 0, n > 2) and another is the high temperature or high chemical 
potential regime (|-u„| ~ 1). We will study the first regime in section 3 and the 
second one in section 4. 

Before proceeding the analysis of our model, we remark about the validity of 
the 1/ D expansion. In the strict large D limit, our analysis will be valid for any 
temperature and chemical potentials. However, in a large but finite D case, some 
problems about the 1/ D expansion will arise in a very low temperature regime and 
a very high chemical potential regime. In the very low temperature regime, the 
dimensionless effective coupling A/T^ will be large. Thus the contributions of higher 
loops will be large and the expansion will not work. In the very high chemical 
potential regime, as we will see later, light mass modes will appear and similarly the 
contribution of the higher loops will become large. We will discuss the details of 
these issues in section 3.2 and 4.2. 

We summarize the regimes which we will consider in this article in Table 1. 

very low temperature regime T/A^/^ < D^'^ 

low temperature and low chemical potential regime l^nl ~ 0, n > 2 

intermediate temperature and chemical potential regime |m„| 7^ 0, n > 2 

high temperature regime |-u„| ~ 1 

high chemical potential regime |-u„| ~ 1 

very high chemical potential regime /i/A^/^ ^ (A;2T/A^/^)^/^ 

Table 1: Various regimes in the analysis of the one dimensional gauge theory. T/A^'^ and 
fi/X^'"^ are dimensionless temperature and chemical potential. We will discuss our model 
in the first three regimes in section 3 and the rest in section 4. 



3. The phase structure of the low temperature and low chem- 
ical potential regime 

3.1 The phase structure in the leading 1/D expansion 

In this subsection, we evaluate the phase structure in the low temperature and low 
chemical potential regime, by analyzing the effective action (2.16). Note that this 
effective action is the leading order of the 1/D expansion. We will consider next 



10 



order corrections in subsection 3.2 but, as we have mentioned, they do not change 
the nature of the phase structure. The obtained phase structure is summarized in 
Figure 1. 



gapped 
gapped? 



non-uiiiforirn 



unknown 




unknown 



Figure 1: Phase diagram of the one dimensional gauge theory in fi — T space from the 1/D 
expansion. Three (uniform, non-uniform and gapped) phases exist and the orders of the 
phase transitions between them are second and third. In the shaded regions, it is difficult 
to analyze the model through the 1/D expansion. In the horizontal shaded region (very 
high chemical potential region fi/X^'^ ^ (^22^/-^^)^)) the expansion does not converge 
because of the existence of the light mass modes. In the inclined shaded region (very low 
temperature region T/A^'^ < D~^), the expansion is not valid, since the effective coupling 
A/T^ becomes too strong. The analysis in the vertically shaded region is also difficult, 
since the Wilson loop operators are highly interacting (intermediate region). However we 
can guess that the vertically shaded region will be gapped phase. See Table 1 also. 



If both of the temperature and chemical potential are sufficiently low such that 
all the coefficients of \un\'^ in (2.16) are positive, the stable solution is given by 
\un\ = for all n. Thus the contributions of |-u„| are small in this regime and it is 
indeed enough to keep only |-ui| in the saddle point equation (2.17) to analyze the 
thermodynamics in this regime. 



T 



l + 2e" 



-/3A 



h + k2 



-^fi 



+ e 



/9m 



Fi 



(3.i; 



Since e ^^, e ^'-^ ^^ and e '^(^+^) will be small in this regime, we can solve this 
equation approximately and obtain the condensation as 



^ = 1 + -e-f^'"'' Ih + h 



-Pix 



+ e 



I3ti 



m 



(3.2) 
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Then, by putting this solution into (2.16), we obtain an effective action for the Wilson 
loops 



DAr2 



where the coefficient a{P,fi) and b{P,fi) are given by 



n=2 



(3.3) 



aiP,^i) = --e-^'''''{h + k2 



b{/3,fi) 



/3A1/3 



-2/3A1/3 



ki + k2 



-/^M _|_ g/^M 



g-/3M ^ g/3M 



(3.4) 
(3.5) 



3 V V 2 

Note that b{l3,^) is always positive. In this case, three phases will appear as we will 
show soon through the argument of the Landau-Ginzburg type analysis in [42, 43]. 
(If b < 0, two phases will appear instead of them.) Here the order parameters of 
these phases will be the values of the Wilson loop operator «„ or equivalently the 
eigenvalue density of the gauge field Aq, which is defined as 

a - a.) = — I 1 + T M„,e-^"'^" 1 . (3.6) 



1 ^ 




Pf 



unifoim 



Pf 



non-uniform 



gapped 




5 

Figure 2: Plots of eigenvalue density function p{a). Three configurations of p{a) charac- 
terize the three phases in the /U — T phase diagram. 

Now we investigate the three phases. If the temperature and chemical potential 
are both low, a(/3,/i) > is satisfied. Then, from the effective action (3.3), the 
stable configuration is |m„| = for all n as we have mentioned. This phase is called 
uniform phase, since the eigenvalue density p{a) is constant and thus it is uniform 
with respect to a. (See Figure 2.) This phase is an analogue of the confinement 
phase in the higher dimensional gauge theory, since the expectation values of the 
temporal Wilson loops vanish. 

As the temperature or chemical potential increases, a{P,fi) becomes 0. On the 
curve a(/3, /i) = 0, stiU |-u„| = is stable since b{/3, /i) is positive. However if a(/3, fj,) < 
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0, the solution |-ui| =0 becomes unstable and a stable solution, which is given by 



|-ui| = a/— a/26 and |-u„| = for n > 2, appears. The configuration of the eigenvalue 
density p{a) becomes non-uniform as shown in Figure 2 and this phase is called 
non-uniform phase. Thus a phase transition happens on a(/3, fi) = 0. This line is the 
first phase transition line in /i — T phase space. (See Figure 1.) It is easy to show 
that this transition is second order by evaluating the free energy [32, 42]. 

As the temperature or chemical potential increases further, |-ui| achieves 1/2. 
Then a gap appears in the eigenvalue density p. (We can choose ui real by a gauge 
fixing, and then the gap arises at a; = ±7r//3. ) This is a Gross- Witten-Wadia type 
third order phase transition [44, 45]. In this phase, all the Wilson loop operators 
become non-zero. This phase is called gapped phase and is an analogue of the 
deconfinement phase in the higher dimensional gauge theory. The curve |-ui| = 



a/— a/26 = 1/2 gives the second phase transition line in Figure 1. 

We have found the two phase transition lines between the three phases. Let us 
analyze the details of these two curves. First we evaluate the curve described by 
a(/3,/i) = 0. For small /x, we can solve this equation and obtain T as a function of 



Tciip) _ 1 ^2 / /^ V^... C3 7A 



f^r 



Ai/3 logD 2 VAi/3 

Here we can see that the existence of the chemical potential reduces the critical 
temperature. It means the chemical potential enhances the non-uniform phase. Note 
that, at /i = 0, this curve is coincident with the result in [32]. 

For finite p, the e^'^ term in (3.4) is dominant and the curve is described as 



Tclif^) _ 1 (, ,,Ml/3 



Thus as p/\^/^ approaches 1, the critical temperature goes to 0, and pd/X^^^ = 1 
seems a critical chemical potential at T = 0. However the 1/D expansion will not be 
valid in such a very low temperature regime and this result is not reliable. Indeed 



®Since the mass dimension of A is 3 in our model, T/A^/"^ and n/X^^^ can be regarded as 
dimensionless temperature and chemical potential. 

^The critical temperature (3.7) goes to zero if we take D -^ oo. Similarly the second critical 
temperature (3.9) will also go to zero in this limit. (In this limit, even though the critical tempera- 
tures become very low, the 1/D corrections are suppressed further and these results are exact.) As 
a result, only the gapped phase appears in D = oo. See Figure 4. 
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some strange things will happen around this point when we consider 1/ D corrections 
in subsection 3.2. 

Next we evaluate the second phase transition line |-ui| = a/— a/26 = 1/2. For 
small /i, the curve is given by 



Tc2{^J) Tclifl) I , , '2 I T^^2 f f^ 



2^ 



,1 + — \l + U\ogDy{J—] +■■■. (3.9) 

AV3 \i/3 y 3D y ^^ ^ ' Vav3 ^ \ \^ ^ > 

For finite /i/A^^^ (< 1) , the curve behaves 



Ai/3 Ai/3 V ^Di-^/\i/3 

However this equation is not valid around /i/A^/^ ~ 1. In order to investigate this 
region, it is convenient to evaluate the saddle point equation for ui, which is derived 
from (2.16), as 



■l _0^ / , e-f'^ + e^^' 



Mi = 0. (3.11) 



Thus in case mi 7^ 0, A has to satisfy 



3A / . . . o-'^^' + e'^^'\ 1 



e-/^- ^A;i + k2 j = -. (3.12) 

Since (3 will be large on the curve near /i/A^^^ ~ 1, we can approximately solve this 
equation as 

A = fi+^logD + ---. (3.13) 

By putting it into the saddle point equation (3.1), we obtain 

The positivity of |-ui| requires that this solution is valid only if A^/\^^^ > 1. Now we 
can derive the second phase transition line around /x/A^/^ ~ 1 by putting |mi| = 1/2 
in this equation, 

^ = l-/x/A^/-Vi + ---. (3.15) 

AV3 ' 6D 

Although this equation predicts a critical value of the chemical potential fj,c2/^^^^ = 
1 + 1/6D at T = 0, the 1/D expansion will not work there. 
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Finally let us confirm that the relation A > fi, which we have assumed, is always 
satisfied in the uniform and non-uniform phase. In the uniform phase, A/A^^^ = 1 
from (3.2). Since the uniform phase exists up to /id/A^/'^ = 1, the relation A > /i is 
satisfied. In the non-uniform phase, the equation (3.12) and (3.13) show A > /x. 

A problem is the case A = fx, which arises on a line T = 0, fi/\^'^ > 1. It causes 
zero modes of the adjoint scalar <l>^ in (2.13). In addition, the 1/D expansion itself 
is not valid around this very low temperature regime. Therefore further analysis is 
necessary but we do not consider it in this article. 

The analysis in the gapped phase is difficult since all the Wilson loop operators 
are excited and interacting each other through a constraint p{a) > 0. (The vertical 
shaded region in Figure 1.) We can perturbatively analyze it just above the curve 



a/— a/26 = 1/2 by assuming that \un\ {n > 2) are small [32, 42]. Thus it is com- 
plicated to show the relation A > /i in general. On the other hand, if temperature 
or chemical potential is enough high, |-u„| ~ 1 is satisfied and we can evaluate the 
contribution of Aq perturbatively. There, analysis is possible as we will see in section 
4 and we can guess the stability problem of the gapped phase through these results. 

3.2 1/D corrections and problems in very low temperature regime 

In this subsection, we evaluate the sub leading 1/D corrections to the effective action 
(2.16) in the low temperature and low chemical potential regime. Then we will show 
that the 1/D expansion is not valid in a very low temperature regime. After that we 
argue how the 1/D corrections modify the phase structure derived in the previous 
section. 

We show the calculation of the 1/D corrections in appendix B and, by using it, 
we obtain the relevant terms of the effective action as 

5(A, {un})/iDN^) = Co + C2IM1P + C^u^l^ + ■■■ + Oil/D^), (3.16) 

where 

^ _ /3A4 PA 



(3.17) 
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C2 =j: - X ( /ci + k2 j + —X ( ki + k2 



A 



A 



A 

4A^ 



4A3 



4A3 



-4 



A 



A 



2-1 



1 + 



4A'^ 



4A3 



4A3 



(3.18) 



Ca =-r-^x'^ I ki + k2 



2D 



A 



4A3 



1(1 A 

2 ^4A3 



(2 + /3A) 



4A» 



/3A ,/ y-y-^\-/ A 



2D 



x<^/3A 



4A3 



4Aa 




Here x = e"'^^ and y = e~^^. Note that higher order terms of x are irrelevant in low 
temperature, since the transitions happen around x ~ 1/D in /i = case and the 
critical temperatures will decrease as /x increases. 

Now we derive the effective action for the Wilson loop operators as in the previous 
section. By solving the saddle point equation for A, we obtain the condensation 



- — = 1 H H — X [ ki + ^2 ] \ui\ + 

AV3 D I 30 32/ 3 V 2 ' ' ' 



(3.20) 



where x = e ^^ . By substituting this solution into (3.16), we obtain the effective 
action for the Wilson loops as 



S/{DN^) = fX^'ho + a>i|2 + 6>i|4 + . . . , 



(3.21) 
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with 



y=.Mt, + fe!^±JLiV ^ + ^/'aV3/2M 2V5V3181 391V5' 



2 J 3 D\\300 9/ 2400 1800 

_ AV^,. f fel^ZiLiy Uv^p^ + ^ + 1 V (3.24) 

£) V 2 y \^ ^400 160 32^ ^ ^ 

From these expressions, we immediately find that the 1/ D expansion in T/A^^^ < 
1/ D regime is problematic, since several 1/D corrections involve /3A^/^ factors. Thus 
if /3A^/^ ~ D (T/A^/^ ~ 1/-D), these corrections become the same order to the leading 
terms. Similar terms may arise in higher 1/ D corrections also and the 1/ D expansion 
will not be reliable in such a very low temperature regime. Thus the 1/ D expansion 
would be valid until T/A^'^ ~ D~^ ^ where 7 is a positive constant^" ^^. This regime 
is depicted as the inclined shaded region in Figure 1. 

Now we evaluate the effective action (3.21) and argue the low temperature phase 
structure involving the 1/ D corrections. We can show that h' is positive at a' = 0. 
Thus the argument in the previous section is still valid and the phase structure does 



not change. The curve a' = and a/— a'/26' = give the two phase transition lines 
in the fi — T plane. However we can see that the curve a' = ends at fi/\^^^ = 1 on 
T = again. This result is the same to the leading result in (3.8). This strange fact 
also indicates that the 1/D expansion is invalid in the very low temperature. 

3.3 General chemical potential 

Until now, we have studied the phase structure only for the simple chemical poten- 
tial (2.3). In this subsection, we take each fij arbitrary value and consider general 
chemical potentials. 



"'^'^7 will be less than 1 through the arguments in this section. In order to determine 7 precisely, 
we have to evaluate the higher order corrections of the 1/D expansion and it has not been done. 

^^In /i = case, the very low temperature regime is in the uniform phase and physical quantities 
do not depend on T. Thus this problem was not observed in [32]. However we cannot rule out a 
slight possibility that new phases appear in this regime. Fortunately, numerical analyses show that 
such new phases do not arise [31, 38, 39, 40, 41]. 
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Before evaluating the general chemical potentials, we discuss a small k2 ~ 1/D 
case, in which the contribution of the chemical potential will be comparable to the 
1/D corrections derived in the previous subsection. Even in this case, the phase 
structure in the low temperature and low chemical potential regime is similar. On 
the curve a' = 0, if /x sufficiently closes to A^^^, k2e^^^'^ ^^^ will be dominant in (3.23) 



even if k2 is small. Thus the /i dependence is still large there. The curve ^ —a! j^h' = 
also depends on /i strongly in a certain regime. Therefore the qualitative nature of 
the phase structure is not modified. 

Now we consider the general chemical potentials. In this case, the one-loop 
contributions from the complex adjoint scalars are modified, and the second line of 
the effective action (2.16) becomes 



^-E 



ra=l 



e 



-nfS 



^ I D-2 [D/2\ ^ 2 ^Y^^ e-"^^^ + e"'^^^ 

^ 2 

1=1 



D D 



\u P 



n 



(3.25) 



During the derivation of this potential, we have assumed A > fij for all fij. 

Here we consider the phase structure. As we have argued in the case of k2 ~ 1/-D, 
even if each contribution of the chemical potentials appears with the 1/D factor, the 
largest chemical potential will be dominant in some regimes to fix the phase structure. 
Therefore we will obtain similar phase structure. Besides we can confirm that the 
condensation satisfies A(/3, {/x/}) > fii in the uniform and non-uniform phase. 

In principle, there is a possibility that a different phase structure appears, since 
we have deformed the potential through the chemical potentials. However it does 
not happen. Especially, in the leading order of the 1/D expansion, we can prove 
that the phase structure is determined by the |-u„| independent terms of the effective 
action (2.16). Since the |-u„| independent terms of the potential (3.25) are the same 
to the previous ones in (2.16), the phase structure does not change. We show it in 
appendix A. 

As we have mentioned in the introduction, it was supposed that the study of 
the low temperature thermodynamics of the gauge theories with the finite chemical 
potential is difficult because of the perturbative instability of the massless adjoint 
scalars. However, in our study, the condensation A(/3,/i) of the adjoint scalars pro- 
tects our model from the instability and we can investigate the low temperature phase 
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structure. In the next section, we will explore the properties of the condensation in 
different regimes. 

4. Condensation in high temperature and high chemical po- 
tential regime 

4.1 Condensation in the large D limit 

We investigate the natures of our model in the high temperature regime and the 
high chemical potential regime. In these regimes, the phase will be the gapped 
(deconfinement) phase and we can use an approximation |-u„| ~ 1 [13, 32, 42]. Then 
we can use a perturbative analysis in Aq around Aq = 0. 

In this subsection, we consider the simple chemical potential (2.3) and we will 
show that a unique condensation A(/3,/i) > fi exists for an arbitrary value of yU in 
the large D limit. This conclusion will be modified in the finite D case as we will see 
in the next subsection. 

In the large D limit, we can ignore fluctuation of Aq and it is enough to evaluate 
the saddle point equation (2.17) only. In |-u„| = 1 case, this equation becomes 



^=l + 2Ve-"'^^ 



n=l 



ki + K2 



2 






The left hand side of this equation is simply increasing from to infinity with respect 
to A. On the other hand, the right hand side is simply decreasing from infinity to 
in A > /i region. As a result, this equation has an unique solution A(/3, /i) in A > /i 
region. (See Figure 3. Similar behaviour can be seen in rf = 2 and 3 dimensional 
gauge theories also.) There is another solution in < A < yU. However this solution is 
unphysical since we have assumed A > /i when we derived the saddle point equation 
(4.1). 

Now we show solutions of (4.1) in several cases. In ki = case, if /3 is sufficiently 
small, we obtain. 



(/3A^/3 < 1). (4.2) 
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In /ci 7^ case, if /3 and fi are small, we obtain 

^ = (j) ^ + ^/j)+--- (/^^/3/A « 1, /3A^/^ « 1). (4.3) 

These two solutions are valid in very high temperature. If A is close to fi, we obtain 

A =/i + ^ + . . . (^/aV3 » (^k,T/~Xy'y/^). (4.4) 

This solution is valid in very high chemical potential regime. Note that this conden- 
sation induces a light effective mass A^ — /i^ ~ 2A;2A//3/i^ for ^^ in (2.13). In the 
next subsection, we will see that such light mass modes cause a divergence in 1/D 
corrections and the 1/D expansion does not work there. Thus the arguments in the 
very high chemical potential regime will be valid only in the D — )■ oo case. 

4.2 1/D correction in very high chemical potential regime 

In this subsection, we evaluate the 1/D corrections in the very high chemical potential 
regime (/i/A^^^ ^ {k2T/\^^^Y^^), in which the equation (4.4) is satisfied, and argue 
the validity of the 1/D expansion. 

In the very high chemical potential regime, since /3(A — fi) will be small and «„ 
will be close to 1, the zero-mode of the Matsubara frequencies of <l>^ will be dominant. 
Therefore the relevant 1/D corrections arise from loops of <I>^ and bab and the path 
integral of Aq. 

First we evaluate the 1/D correction from Aq. In the very high chemical potential 
regime, the dominant A dependent terms in (2.15) can be evaluated as 

hog {{a^ - a,)' + 2 {A' + fi') {a, - a,f + {A' - fi'YY (4.5) 

Then the effective action for Aq is given by 



^log{{a,-a,)' + 2{A' + f,'){a,-a,f+{A'-f,'Y]-^\ogi 



(4.6) 
where the last term is derived from (2.11) by assuming that ttj are small. Now we 
assume A^ — /i^ ^ {aj — UiY. Then we can expand the first log term and obtain^^ 

N'D log(A^ - ^^^) + 2ND 2,_% E "' " E ^ l°g(". " "^)'- (4-7) 

^ ^ ^ i=i i,j 



^^li we start an assumption A^ — ij,^ -^ {aj — ai)"^ in (4.6), the attractive force between a^ will 
be quite strong and this assumption will not be satisfied. 
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Note that this action is just a gaussian in Aq, 

N'D log(A^ - /.2) + 2\hi^±jQ-TTAl/g'. (4i 



Since the coefficient of A^ will be enough large in jj,/\^^^ ^ {k2T/\^^^Y^* regime, ai 
will be strongly trapped around ttj = 0. Thus the assumption A^ — jjp' ^ {aj — aiY 
will be satisfied. Then the gaussian integral of Aq gives a 1/ D correction and the 
effective action for A in the very high chemical potential regime will become 

S,,(^)IDN'' ^ - i; + § .og(A^ - , V ^ log (g^) 

+ I O I — J from matter loops )+■■■. (4.9) 

Here the first and the second terms are from (2.16) with an approximation /3(A — 
/i) ^ 1. The third term is from the Aq integral. Thus the 1/ D correction from Aq 
is qualitatively not important for large D case. However, if D = 1, this correction 
cancels the second term. As a result, the arguments in the previous section will not 
be valid and (4.4) will not be satisfied. Thus the assumption A ~ /i in /i/A^^'^ ^ 
{k2T /X^^^y/'^ is not ensured. Therefore a different analysis is necessary in D = \ 
case and we will discuss it in the next section. 

Now we evaluate the sub leading 1/ D corrections from the matter loops. In the 
large D case, we have confirmed that Aq is sufficiently small and we can ignore it in 
the loop calculation at this order. Then the dominant contributions of the matter 
loops are from the zero modes of the Matsubara frequencies of ^^ and hab- Therefore 
we can evaluate them by using a zero dimensional reduced model 

Sod = -^kbM^b'Jcd + Yl ((^' - /^')*i'*a + g^bab^i'H) ■ (4.10) 

1=1 

Here bat and $^ are the zero-modes of the one dimensional fields. Then we can 
calculate the 1/D corrections of the effective action (4.9) as in appendix B of [32], 

1 A 1 / k2\T \"'\ 

From this expression, we notice that if (A^ — yU^)^ < k2XT, the last sum does not 
converge. It means that the 1/D expansion does not work in this regime. From (4.4), 



21 




it will happen 

fi* > 4k2~XT. (4.12) 

Note that this estimate is crude, since, if /i is not sufficiently larger than (A;2AT)^/^, 
the reduced model analysis (4.10) is not valid. Especially in the uniform and the 
non-uniform phases, the contribution of the gauge field is relevant and the 1/D 
expansion still works as we have discussed in section 3.2 even if (4.12) is satisfied. 
By considering it, we conclude that the 1/D expansion is not valid in the horizontal 
shaded region in Figure 1 schematically. Although the 1/D expansion does not work 
in this regime, we cannot conclude that the system is unstable there. The divergence 
of the 1/D correction in (4.11) arises from the loops of the light mass modes of ^^ 
and it is not clear whether it indicates an instability of the system or not. 

1/D correction in finite chemical potential Here we consider the 1/D cor- 
rections in the finite chemical potential regime fi/X^'^ < {k2T /\^''^Y'^. In a finite 
temperature, the calculation for the corrections will be complicated but, if temper- 
ature is enough high, a zero dimensional analysis similar to (4.10) is possible. Then 
it is not difficult to show that the 1/D corrections converge if (A^ — /i^)^ > k2XoT 
is satisfied. Therefore we can guess that the 1/D expansion is valid in the finite 
temperature case also. 

By using the 1/D expansion method, we have revealed that, if the chemical 
potential is not very high, our model in the high temperature regime and the low 
temperature regime is stable. The stability in both of the regimes will support the 
stability of the unknown regime in the gapped phase. (The vertically shaded region 
in Figure 1.) However we have investigated only one condensate vacuum and there is 
a possibility that this vacuum is just a local minimum and more stable vacua exist ^^. 
Another possibility is that the model is unbounded below in the finite chemical 
potential. Thus we conclude that our model is at least meta-stable if the chemical 
potential is not very high. 

^■^The appearance of other phases might be natural If fi is large, (Tr|$^p) ^ (TrF*^) wih be 
satisfied since the chemical potential makes the effective masses of $^ light. Then the eigenvalue 
distribution of the adjoint scalars will be pancake like, and, intuitionally, a doughnut like distribution 
may be favoured. If such a transition happens, it may correspond to the Meyers-Perry black 
hole/black ring transition in higher dimensional gravity [11]. Besides, several intermediate deformed 
Meyers-Perry black hole solutions also exist in gravity [46]. However, we have not found such new 
phases in our model and it is interesting to investigate them further. 
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4.3 Condensation in D = 1 case 

As we have seen in equation (4.9), tlie D = 1 case is special. In this subsection, 
we show that a consistent condensation will happen in D = 1 case also by using 
a different analysis up to 1/D order. However, this condensation will give rise to 
a small effective mass A^ — /i^ in a very high chemical potential and it will cause 
a divergence in the next order of the expansion. This is similar to the behaviour 
of the D > 1 case and the 1/D expansion will not work in the very high chemical 
potential regime. In addition, we will show that our model has complex fuzzy sphere 
like saddle points in Z^ = 1 case, although their physical interpretations have not 
been understood. 

For simplicity, we consider a sufficient high temperature regime {(3\^^^ <^ 1) 
only. Then the zero-modes of the Matsubara frequencies will be dominant and the 
model reduces to a zero dimensional model, 

S =Tr {-2i2igo<^^ [Ao, $] - /i^$^$) 

+ 5^ gl (^$1$, + ^-Y:Y,^ + ^AoaAo,^ M,,,,, (^$1$, + ir/rj + ^-Ao.A^,^ , 

(4.13) 

where we have scaled matrices appropriately and the coupling is defined as g^ = g^T. 
By employing an auxiliary matrix Bab = ^"^^ab + Qobab where hab satisfies haa = 0, 
this action becomes. 



S =Tr ( -2^,g,W [A,, X] + ^-^ {X' + W) + ^^^ + ^ ^^'" 



i=l 

r2 A4 



DN /\ 1 a 

. -habM-,\^Kd + ^hab {Y^Yi + XaXb + WaWb + A^aA^b) ■ (4.14) 



Here we have used $ = {X + iW)/\/2 and Aq = g^ND. We evaluate this action up to 
the second order of the 1/D expansion. In this case, we can ignore the interactions 
between Aq^X,W and hab in the last term. By integrating out Y^ and hab through a 
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technique in [32], we obtain 
S=Ti(-2figoW[Ao,X] + 



A^ - //2 A^ 



+ DN^ 



1 



8An 4 



' A4 



1 /Mo 

2 I A^ 



+ 7^ log 1 + 




(4.15) 



Here ki = {D - 2)/D. 

Now we evaluate the path integral of X, W and Aq and derive an effective action 
for A. A formula for the following three matrix model is available [47]: 

5 = Tr |aMi[M2, M3] + ^ (M^ + M^) + ^mA , 
DM1DM2DMS exp (-5) 



Z 



=Ca ^ / (imi ■ ■ ■ dm]\f TT 



m,- — m 



rrii — rrij + 



in 



(4.16) 



where a, 6, c are constants and Am = N/gj^ = dp' /2a^ . C is an irrelevant factor. We 
have ignored 1/N corrections. Here the free energy is given by 

1. . 1 



Fq{\m) ^ -| log^^ + -git + 



-)■ 



3(127r)2/3 _2/3 



40 



-9m 



0{g 



-5/3n 
M ) 



(Qm -^ 0), 
{gM -^ 00). 



In our case, from (4.15), gM becomes 



9m 



ifi'^Xo 



(4.17) 
(4.18) 

A2(A2-/i2)2/}- ^^-^^^ 

If /i is small, gM will be also small. Then (4.17) gives us standard log terms plus 
0{1/D'^) corrections. Thus, in the small fi case, a consistent condensation A > /i 
will happen as usual. On the other hand, if /i is large and A^ — /i^ is small, then g^^ 
will be large. In this case, we obtain the effective action for A as 
S(A) A"^ '" Q/'1o^^2/3 /a2,'a2 ,,2^2n^V3 

°'\o 



DN^ 



Un 4 °^ 40D V 8/i2Ao 




2 log I 1 




+ ■■ 



(4.20) 
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From this action, we can derive a saddle point equation for A and we can obtain a 
unique saddle point A(/i) which satisfies A(yu) — /i ^ for any fi. 

Therefore, up to this order of the 1/D expansion, we obtain the consistent con- 
densation and the system is stable even in D = 1 case. We consider the high 
temperature approximation in this section but similar analysis will be possible in 
a finite temperature also. However, since A^ — yU^ will be small in the very high 
chemical potential regime, the next order terms will diverge and the 1/D expansion 
will not be valid as in the large D case. (The appearance of the fractional power of 
D in the third term of (4.20) also implies a problem of the 1/D expansion in this 
regime.) 

Fuzzy solutions? Now we discuss possible saddle points of the zero dimensional 
action (4.15). Since the cubic interaction in (4.15) can be regarded as a CS like term, 
the action has a complex fuzzy sphere like saddle point [48], 

X = -I Ji, W = -I J2, Ao = -t^ J3, (4.21) 

2/i5f 2(ig 2^g 

where Jj are the generators of the A^ dimensional irreducible representation of SU(2), 
which satisfy [Ji, Jj] = ieijkJk- By replacing Jj with reducible representations, we 
can obtain many saddle points. However, these fuzzy sphere like solutions are not 
hermitian^^ and physical interpretation is unclear. 

In addition to these complex saddle points, it might be possible to find different 
fuzzy solutions in (4.13) or (4.15) as in the studies in [34, 48, 49, 50, 51]. 

General chemical potentials Now we consider the general chemical potentials 
as in section 3.3. If the values of several chemical potentials are the same and larger 
than the others, the analysis in the previous section is valid. If only one chemical 
potential is very large, we can approximately apply the analysis in this section by 
ignoring other chemical potentials. Thus the consistent condensation A > yU/ will 
always happen up to the 1/D order. 



-'^^At these saddle points, A < /i may not be forbidden. (However it will be unstable.) Then X 
and W can be hcrmitian but ^0 is still not. 
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5. High temperature condensation in higher dimensional gauge 
theory 

In this section, we consider the generahzation of our argument about the condensa- 
tion to d dimensional gauge theory, 



S 



/3 



dt / d 



jd-l. 



X 



D 



^9l ''' 



1=1 



I,J,'Lj ^ ^ ^ 




(5.1) 



Here gd is the gauge couphng. We consider the simple chemical potential (2.3). We 
will show that the condensation of the adjoint scalars can happen at least in the high 
temperature regime in the leading order of the 1/ D expansion. 

First we employ an auxiliary field Bab as in the d = 1 case. Then we assume 
that this field condenses as Bab = ^"^^ab-, where A does not depend on the time 
and position. We also assume that this condensation satisfies A > /i. Under these 
assumptions, we can exactly derive a saddle point equation to determine the con- 
densation A in the large D limit. 

In the large D limit, we can ignore the contribution of the spatial components of 
the gauge field Ai by regarding D ^ d—1. The interactions between hab and ^^ and 
y* are also suppressed. (Here hab is defined as in (2.9) and satisfies J dtd'^~^xbaa = 0.) 
In the high temperature regime, we can also ignore Aq. Then, we can integrate out 
^^ and y* and obtain the effective action for A^ as in section 2.2. From this effective 
action, the saddle point equation is obtained as 



A^ 



E 
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d' 



d-i 
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(27r: 



id-l 
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2ki 



( 2Trn 



+ p2 + A2 



k. 



fe+2/i) +p2 + A2 f^-^/i) +p2 + A2 



(5.2) 



Here \d = g^ND is the d dimensional 'tHooft like coupling. Generally solving this 
equation is still complicated and we evaluate it by taking a high temperature limit. 
Here the zero modes of the Matusbara frequencies are dominant and we can ignore 
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the non-zero modes. Then the equation is simphfied as, 



(27r 



v(i-l 



2K 



d-2 



A 



5-d 



\,T 



A2 

A2 



=A;i/rf(AVA2) + A;2/rf((A2-/,2)/A2 



(5.3) 



where 



fd{x) = A 



A „d-2 

g P^ + XA^ 



(5.4) 



Here Ki_2 is the volume of the d — 2 dimensional unit sphere and A is a momentum 
cut oflF. T appears only through A^T in this equation. Thus we can regard X^T as 
an effective coupling. This equation determines A in terms of the effective coupling 
XdT, fi and the cut off A. We tune A dependence of X^T such that the condensation 
A is fixed for a particular physical value at a certain chemical potential /x. 





Figure 3: Solution of (5.3). In d < 3 case (the left plot), the curve represents the schematic 
behaviour of the right hand side of (5.3) in A^ > ^^ region. It diverges at A^ = fi'^. The 
straight line is the left hand side of (5.3) and the crossing point gives the solution of (5.3). 
Note that a unique solution exists for any value of /U. In d > 4 case (the right plot), the 
curve does not diverge at A^ = A*^- As a result the curve cannot cross the line for small 
Xd (or large ^) and the solution does not exist in this case. 

Now we evaluate the equation (5.3). We show the explicit expressions for fd(x) 
in appendix C and we can derive the condensation by solving it. Instead of doing 
it, here we argue the condition for the existence of the condensation in A > /i 
region through the qualitative properties of fd{x). We can show that fd{x) is simply 
decreasing in a; > and behaves as 



-^ oo X -)■ +0 {d<3) 

fdix){ ^l/{d-3) x^+0 (rf>4). 
-)> X -)■ +00 



(5.5) 



Thus the right hand side of (5.3) diverges at A = /i in the d < 3 case but does not 
in the d > 4 case. Here the consistent solution of (5.3) is given by a crossing point 
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in Figure 3. Therefore the equation (5.3) always has a unique consistent solution in 
d < 3. On the other hand, in li > 4 case, one solution exists only if the following 
equation is satisfied. 



\dT ^ (27r 



id-i 



A5-d 



fl 



(5.6) 



2\/rf_2 VA; A;i/,(/iVA2) + A;2/(rf-3)- 

We have obtained this condition by evaluating (5.3) at A = /i. Thus the consistent 
condensation does not happen if the effective coupling is weak. Equivalently we can 
say that, for a given coupling, a critical chemical potential Hc exists, which saturates 
(5.6), and the condensate happens only ii fi < n^- We summarize this result in Figure 
4. 



T 



d=l 



T 



gapped 



rf = 2,3 



T 



gapped 



unknown 



l-tc d>4 
' '^ ^ ' / unknown 



unknown 



Figure 4: Schematic phase diagrams of the d dimensional large N gauge theories in ^u — T 
space at D = oo. In the D = oo case, the problems from the 1/D corrections do not arise. 
In d = 1 case, since the two critical temperatures become at D = oo, only the gapped 
phase appears. (See (3.7) and (3.9).) In d > 2 case, we evaluated only the high temperature 
regime. In d > 4 case, there is a critical chemical potential and the condensation does not 
happen beyond it. The nature of this regime has not been understood. 

One important unsolved issue is the nature oi fi > fic region, in which the 
condensation does not happen. One possibility is that the system is destabilized 
by the chemical potential. Another possibility is that the system is still stable but 
described by a highly interacting theory. 

Although the above results are exact in the large D limit, as we have observed 
in section 4.2, the 1/D corrections are important in the very high chemical potential 
regime in finite D case. Thus our results may be valid only in the D — )> oo case in 
this regime. 

Even if we start the supersymmetric gauge theory^^, the behaviour will be the 
same in the high temperature regime, since all the fermions are decoupled if the 



"'^^How to take large D limit in the supersymmetric gauge theory is difficult problem, since the 
number of the bosons and fermions grow at different rates as D grows large. By taking the high 
temperature limit and ignoring fermions, we can ignore this problem. 
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temperature is sufficiently high. Therefore the analysis in this section may be valid 
in the supersymmetric gauge theory also. 



6. Conclusions and discussions 

In this article, we investigated the thermodynamics of the large A^ gauge theories 
with the chemical potentials. Because of the condensation of the adjoint scalars in 
the large D, we can analyze the effect of the finite chemical potential even in the per- 
turbatively massless gauge theories. It is an important step towards understanding 
the phase structure of large A^ gauge theories. However the light mass modes of the 
adjoint scalars in the very high chemical potential regime cause the 1/ D expansion 
to diverge. Understanding of the nature of this divergence is important to figure out 
whether the system is stable or not in this regime for finite D. In addition, if we can 
understand the 1/ D corrections in the higher dimensional gauge theories, it may be 
possible to apply our analysis to D brane theories at a high temperature. It would 
then be interesting to compare our results in the weak coupling regime to the strong 
coupling results predicted by the dual gravity analysis [3]. 

We also studied the phase structure of the one-dimensional gauge theory as in 
Figure 1. As we have mentioned in the introduction, this theory is related to Dl 
branes on a small circle and it is interesting to evaluate this system by using the dual 
gravity and compare each other. 

However we have not completed the understanding of the phase structure of our 
model. There is a possibility that other phases appear. If such phases are found, 
they may correspond to non-spherical black objects like a black ring or hairy black 
holes in the dual gravity. 

In the sufficient high temperature regime, the one-dimensional model reduces to 
the zero-dimensional model (4.13). This model is similar to the bosonic IKKT ma- 
trix model with negative mass and imaginary CS like terms. As we have discussed in 
section 4.3, fuzzy solutions might exist in this model [34, 48, 49, 50, 51]. Therefore, 
exploring this matrix model may be the simplest way to find the new phases. 
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Our phase structure in the bosonic gauge theory is quite different from the re- 
sults in the one-dimensional supersymmetric gauge theory predicted from DO brane 
black hole [1, 3]. There, the low temperature phase transition does not happen and 
the system is always in the gapped (deconfinement) phase. Besides if the chemical 
potential is larger than a critical chemical potential Hc = cT, the system is desta- 
bilized. These differences indicate that the contribution of the fermions is relevant 
in the low temperature regime. One possibility is that, in the supersymmetric the- 
ory, the effective action (2.16) is modified such that the condensation behaves as 
A(T, /i) ~ T in the low temperature regime. Then the potentials for the Wilson 
loop operators may be unstable at |-u„| = even around T = and the gapped 
(deconfinement) phase will be preferred. In addition, the system will be destabilized 
if /i > A ~ T. However, how to take a large D limit in supersymmetric theories has 
not been understood and we cannot show such a mechanism yet. 
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thank Avinash Dhar, Oleg Evnin, Shoichi Kawamoto, Hiroaki Kohyama, Manaven- 
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for useful discussions. 

A. Cominent on phase structure in 1/D expansion 

As we have seen in section 3, the effective action for Un determines the phase structure 
of the one-dimensional large N gauge theory. Especially, the potential for |-ui| is 
relevant in the low temperature and low chemical potential regime [42, 43]. In this 
appendix, we argue how the phase structure is fixed in the 1/D expansion. 

Generally, the relevant terms in the effective potential for the condensation A^ = 
r and Ui will be given as 

S{t,Ui)/DN' = Coir) + C2(r)|Mi|2 + C,{t)\u,\^ + ■■■ . (A.l) 

Thus the saddle point equation for r is 

C'oir) + C',{r)\u,\' + C',{t)\u,\' = 0. (A2) 

Here '"" denotes d/dr. We assume that C2 and C4 are small compare to C'q. Indeed, 
by estimating the 1/D corrections, we will see that Cq, C2 and C4 are order 1, 
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1/ D and 1/ D^ respectively near the critical point. Then we can solve this equation 
perturbatively by expanding r = tq + ri + ■ ■ ■ . First tq is given by a solution of 



Next Ti is given by 



Tl 



C^(ro) = 0. 



1 ^~(/// \ 1 



(A.3) 



(A.4) 



By putting this solution into the effective action (A.l), we obtain an effective action 
for the Wilson loop, 



S{u^)/DN' = C^{t^) + C^{t^)\u^ 



C4(t-o) 



2 C''(ro) 



l«i|' + 



(A.5) 



The sign of the coefficient of \ui\'^ fixes the phase structure though the Landau- 
Ginzburg type argument in [42, 43]. If it is positive, three phases (uniform, non- 
uniform and gapped) will appear and the phase transitions between them are second 
and third order. If it is negative, two phases (uniform and gapped) will appear and 
the phase transition is first order. 

In the 1/ D expansion, since C4 is very small, the sign of Cq{tq) is relevant in 
(A.5). In our case, from (3.17), Cg (tq) is a negative constant and the phase structure 
is fixed as the former case. 

Note that Co(t) can be regarded as an effective action for the condensation A in 
the uniform phase. Thus the discussion in this appendix implies that this effective 
action determines the phase structure. This observation is interesting but physical 
meaning is unclear yet. 

B. Calculation of the 1/D correction 

In this appendix, we show the derivation of the effective action (3.16) involving the 
subleading 1/D corrections. The analysis is similar to the one shown in appendix E 
of [32] and we do not show the details here. 

The chemical potential dependence appears only through the propagator of the 
complex adjoint scalars 



,nu\^J^-^ ^ 



i^(aj-ai)+^l)\\t\\ 



i^imiim 



2A 



-Allil 



.Alltll 



oil3(aj~ai)^-l3{A~iJ.) ^ _ Qil3{aj~ai)^l3{A+fi) 



(B.l) 



ij 
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Here ||t|| denotes ||t + n/3|| = t for < t < /3. By using this propagator, we can 
evaluate the 1/D corrections to the effective action from loop diagrams of $^, Y^ 
and bab- 

li ki = {D = 2D), we do not need to consider Y^ and the {m + l)-loop 
correction to the effective action is calculated as [32] 

n=— oo i,j=l 

where dm is a factor derived from a property of Mat cd as, 



0?! = —1, (^2 = 3, dm = i (m > 3). 



Gj^lf^ is defined as 



^n,ik — o A2 \^n,ik^ik + ^n,ik^ik + Qn,ikOQ,ik) , (B-3) 

where n dependent term P^iki Pnik ^^^ Qn,ik are defined as 

1-1 ^^11 

p — p+ 



n,Jfc • i/3(afc-a,)-2A/3 _ ' «,*fc ^V i/3(afc-a,)+2A/3 _ ' 

7^^ iTTZ 

27rj "- 



and n independent term S*^^, S^ and S'^^j^ are 



^+ =1 + ^ X- (r"«-™«t„ + l/'"^™^'-™) , (B.5) 

m=l 

5,7 =1 + ^ X'" (y'"W™M!.„ + y-'^U.mU'm) , (B.6) 

m=l 



k,m=0 



(B.7) 



where a; = e~^^ and y = e''^^. 

li ki 7^ 0, we have to evaluate the contribution of Y^ also, but the same formula 
(B.2) is still available with simple changes, 

Sq,u — > ki Sq^uI^^q + k2SQ^ii, 
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By using technique in [32], we can evaluate (B.2). The results in /i = case is 
given by (E.30)-(E.32) in [32]. In the non-zero chemical potential case, the modifi- 
cations on (E.30)-(E.32) are summarized as, 

• (E.30) is not modified. 

• (E.31) is modified as x ^ x iki + k2^^ — ) only. 

• (E.32) is modified as 






X 



<^« - ^'" ' 2 + Ap] (k. + k.^±f^)' + r&i4^ . A,) ki 'y -y-'^' 



(2n-2)!! ^^y 1^ 1 ■ - 2 ) ' V(2^-4)!! / V 2 

(B.8) 

Then by summing over n and adding the leading results (2.16), we obtain the effective 
action (3.16). 

C. Expressions for fd{x) 

In this section, we show the expressions for fd{x) up to rf = 5, which are important 
to show the existence of the condensation in the higher dimensional gauge theories. 
fd{x) is given by 



Then we can calculate it as 



f2{x) =^1= arctan ill ^/x) , 

\/x 



/3(x)=ilog(l + l/x), 

fA{x) =1 — v^ arctan {\/y/x) , 

hi.^) =2 ~ 2^1og(l + l/x). 



These satisfy the relation (5.5). 



Oix^ 
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